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core repulsion on a cZ-dimensional lattice. In this model, which includes the asymmetric exclusion 
process, contact processes and other processes, the stochastic variables are particle occupation 
numbers taking values = 0,1. We show that on a 10-paramctcr submanifold the fc-point equal- 
time correlation functions ( n^i ■ • ■ ng^, ) satisfy linear differential-difference equations involving no 
higher correlators. In particular, the average density ( ng ) satisfies an integrable diffusion-type 
equation. These properties are explained in terms of dual processes and various duality relations are 
derived. By defining the time evolution of the stochastic process in terms of a quantum Hamiltonian 
H, the model becomes equivalent to a lattice model in thermal equilibrium in d -|- 1 dimensions. 
We show that the spectrum of H is identical to the spectrum of the quantum Hamiltonian of a 
d-dimensional, anisotropic spin-1/2 Heisenberg model. In one dimension our results hint at some 
new algebraic structure behind the integrability of the system. 
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1 Introduction 



An interesting class of non-equilibrium problems with a rich dynamical behaviour are 
stochastic reaction-diffusion systems [Q, ^. These processes may involve one or several 
species of particles A, B,C, . . . and an inert state equivalent to the absence of any of 
the interacting particles. Examples of such processes are, to name but a few, coagulation 
A + ^ ^ A , pair annihilation A + yl — > |0, ^, |6| or two-species annihilation A + B ^ $ 
0,|8|. Formulated as a lattice model, corresponds to a vacancy on a site of the lattice and 
particles represented by particle occupation numbers may hop in the lattice (^0 — > 0^) 
and take part in the reactions. Such lattice systems are, in general, difficult to treat by 
rigorous means and correspondingly, considering the vast amount of such models, relatively 
few exact results are known. 

Over the past few years the formulation of stochastic processes in terms of quantum 
spin systems has turned to be a convenient tool in the study of non-equilibrium lattice 
problems (see e.g. ||5|, |^, 11, 12, 13, 14] and references therein). A paradigmatic 



example is the representation of symmetric diffusion of hard core particles (known as 
the symmetric exclusion process |15, |l^]) by the spin-1/2 Heisenberg model (see Q for 
a detailed discussion). But also asymmetric hopping ^, 10, 11, multimer processes 



1^] and reaction-diffusion processes |6|, |13|, 14 1 have been similarly represented. By these 
means, standard techniques for quantum spin systems such as spin wave theory ||5|, [T7[ |, 
Bethe ansatz and related algebraic techniques |6|, |l^, 11, |l^, 14, 18 1, global symmetries 
1^, 12 1 and Goldstone broken symmetry arguments Q have given many new results for 
stochastic systems. 

An intriguing feature of these models is that even though these are interacting many- 
particle systems, some of them are known to give rise to closed systems of differential- 
difference equations for time-dependent correlation functions and exact results have been 
obtained. In one possible scenario the time derivative of a fc-point correlation function 
does not involve higher order correlators, or, in other words, one obtains a closed set of 
not more than k coupled, linear, differential-difference equations. Some well-known ex- 
amples of single-species processes where this happens are the symmetric exclusion process 
and symmetric partial exclusion process describing diffusion of particles on a lattice in 
any dimension |9|, [l^, |l^, the asymmetric exclusion process describing driven diffusion 
in one dimension |[l^ or the voter model describing annihilation (death) processes and 
decoagulation |l^ (for details see below). In another scenario, certain subsets of correla- 
tion function decouple from each other. This happens e.g. in the Glauber model |19|, in 
random sequential dimer deposition [^] or in the generalized models studied in Ref. |21]. 

This observation raises a number of questions, the most obvious one being whether 
there is a classification of these processes, i.e., a general criterion on the reaction and 
diffusion rates such that the resulting equations for the correlation functions decouple. 
Clearly, the answer depends on which correlation functions one wishes to study. Here 
we consider density correlation functions which are the ones one is usually interested in. 
Other correlation functions e.g. exponentials of integrated densities [ff2l or particle-string 



correlation functions [12, |2l[ give rise to other necessary and sufficient equations for the 
rates leading to decoupled equations. As we shall show below, the quantum Hamiltonian 
formalism that we use throughout this paper is a convenient tool for the derivation of these 



equations ( |3.8| ) for the rates and it opens the way to a (partial) physical and mathematical 
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understanding of this phenomenon. This is because by expressing the process in terms of 
a quantum Hamiltonian one relates a d-dimensional non-equihbrium problem to a d + 1- 
dimensional non-equilibrium problem into which one may have (and in the cases discussed 
here actually has) some insight. Here we study only stochastic processes of hard-core 
particles with a two-site interaction. The corresponding quantum Hamiltonian turns out 
to be that of a generalized spin-1/2 Heisenberg model. However, our strategy is easily 
generalized to many-species models or to models with interactions involving more than 
only two sites. 

In some of the known cases which are contained in our more general model the observed 
decoupling of the correlation functions can be understood in terms of a dual stochastic 
process p|, H]. In the case of the (self dual) symmetric exclusion process duality relates 
the time-dependent A;-point density correlation function to the process with a /c-particle 
initial state. Because of particle number conservation this is a great simplification and 
many exact results have been obtained in this way. Also for other processes duality may 
be used to derive new results [16| and so one other question we discuss is the existence of 
dual processes to those which satisfy the constraints (3.8) on the rates discussed above. 
The dual processes that we shall obtain involve additional restriction on the rates arising 
from the positivity of the dual rates and conservation of probability in the dual process. 
In any case, the dual process contains only hopping and various annihilation terms, but 
no non-zero particle creation rates. This is another way of understanding the decoupling 
of the correlation functions from higher order correlators. 

One more problem that we shall address, albeit only briefly, is that of the integr ability 
of the system. From the structure of the equations derived in Sec. 3 and from the solution 
of these equations for the one-point function (i.e., the density profile) obtained in Sec. 6 
it becomes apparent that the system is partially integr able in any space dimension. This 
means that some of the equations for the /c-point functions are integrable and therefore 
yield the spectrum of a subspace of the Hamiltonian. In Sec. 5 we show that the spectrum 
of H is identical to that of a spin-1/2 Heisenberg Hamiltonian which, in one dimension, is 
completely integrable by the Bethe ansatz. However, the generalized Hamiltonian H is not 
related to the Heisenberg Hamiltonian by a similarity transformation. The integrability 
of H can not be derived from the usual Baxterization procedure [^] and the algebraic 
structure underlying the generalized model remains an open problem. 

The paper is organized as follows. In Sec. 2 we give for the benefit of the reader not 
familiar with the quantum Hamiltonian formalism a full discussion of its relation to the 
usual description of the stochastic process in terms of a master equation. We also introduce 
various definitions and relations used later. We define the problem on a hypercubic lattice 
in d-dimensions with periodic boundary conditions even though many of the results derived 
later are also valid for other lattices or boundary conditions. When appropriate, this is 
indicated by an additional explicit remark. Otherwise one should always think of the 
periodic hypercubic lattice. In Sec. 3 we derive the equations for the rates such that one 
obtains linear, inhomogeneous differential-difference equations, i.e., equations for the k- 
point correlation functions which are decoupled from higher order correlators. In Sec. 4 
we discuss duality and derive the criteria on the existence of a dual process. This leads 
also to a number of duality relations for the correlation functions. In Sec. 5 we discuss the 
mapping of the stochastic problem to the Heisenberg quantum Hamiltonian and in Sec. 6 
we calculate the time-dependent density profile with an arbitrary initial state. In Sec. 7 
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we summarize the main results and present some open questions. 



2 Stochastic processes in the quantum Hamiltonian formal- 
ism 

We study one-species exclusion processes in d dimensions, i.e., a system of particles on a 
hypercubic lattice with M sites where each site is either empty or occupied by at most 
one particle. The state space of the system is therefore X = {0, 1}^^ and a given state of 
the system may be represented by a configuration n = {ni,n2, ■ ■ ■ , um} where nj = 0, 1 
and 1 < i < M labels the sites of the lattice. An alternative possibility is to give the 
set {xi,X2, ■ ■ ■ ,xiy} of occupied lattice sites. In this notation, the empty set represents 
the empty lattice and 1 < < M is the total number of particles in the configuration. 
Xi = {x\^\x\^\ . . . jxf^^) is a d component object defining the (integer) coordinates of the 
particle in the lattice. When working on finite lattices, we shall label each space coordinate 
by an integer 1 < x^") < L^'^\ for an infinite lattice x^""^ G Z. For later convenience, it is 
useful to introduce also the unit vector in a- direction, e^"^ = (0^^), . . . , l^''), . . . ,0('^)). 

The stochastic dynamics of the system may be defined in terms of a master equation 
for the probability f{n;t) of finding the configuration {ni, n2, . . . , n^/} at time t. We 
shall use quantum Hamiltonian language which has proven to be a useful formalism for 
stochastic processes on lattices. Each state n £ X is represented by a vector | n ) (or 
\ xi, . . . ,xn), with I ) = I ) being the empty state) and the probability distribution is 
mapped to a state vector 

l/(i))= E/fe*)ll^) • (2-1) 

The vectors | n ) together with the transposed vectors ( n | form an orthonormal basis of 
|-^2-j®M ^YiQ time evolution is defined in terms of a linear 'Hamilton' operator H acting 
on this space of dimension 2*^ 

^\fit)) = -H\f{t)) . (2.2) 
A state at time t = to + t is therefore given in terms of an initial state at time to by 

|/(to + r) =e-^^|/(to) . (2.3) 
From ( |2.1D and (2^) and using f{n;t) = (n| f{t) ) the master equation takes the form 



§-J(.n;t) = -{n\H\f{t)) . (2.4) 

Note that 

(s|/(i))= E/fe^) = l (2-5) 

where 

{s\ = Y.i^\ (2-6) 
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which expresses conservation of probabihty. This imphes {s\H = for any stochastic 
process. The right eigenvector (s) of H with eigenvalue £"0 = and normalized according 
to ( ^.5| ) is (are) the steady state (s) of the stochastic process. In general H is not symmetric 
which means that the rate w{ri]'nl) = —{n\H\ri!) with which a configuration n' switches 
to a configuration n is not, in general, equal to its reverse rate w{n';n). As a result the 
stationary distribution(s) S{n) may be highly non-trivial. The real part of all eigenvalues 
of H is larger or equal to zero. 

Average values ( Q ) are calculated as matrix elements of suitably chosen operators Q 
which may be expressed in terms of the usual Pauli matrices ai'^'^ acting on site j. A 
complete set of observables are the occupation numbers = 0,1. Defining projection 
operators on states with a particle on site j of the chain as 

"7 = 5('-''f) = (o 1)- 

one finds that the average density of particles at site j is given by {n-j) = {s \nj\ f{t) ). 
Correlation functions ( n^^ ■ ■ ■ n^^ ), i.e., the probabilities of finding particles on the set of 
sites {xi, . . . ,Xk}^ are computed analogously. Note that in ordinary quantum mechanics 
average values would be taken as matrix elements between normalized eigenstates of H, 
i.e., {Q) = {k\Q\k) whereas here an average value is the quantity {s\Q\ f) where | / ) 
is (in general) not an eigenstate, but a state with real coefficients < f{n;t) < 1 (2.1) in 
the basis spanned by the set {\n)} and normalized such that ( s | / ) = 1. 

For later convenience we also introduce the operators = (a^ it io"H)/2. In our 

j J j 

convention 

/ \ 

1 X 

creates a particle at site j when acting to the right, while 

st=(^ ) (2.9) 

annihilates a particle at site j. Note that 

{s\st = {s\nj and {s\sj = {s\{l - nj) . (2.10) 

Introducing the ladder operator = J2"^^ one may write 

^ j 

{s\ = (0|e^^ . (2.11) 



Using the commutation relations for the Pauli matrices then yields ( |2.10 ) 



Now we are in a position to define the stochastic processes we intend to study by a 
quantum Hamiltonian H. We define 

1 
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with the nearest neighbour reaction matrices 

/ an ai2 ais 014 \ 



(a) 
3 



«21 «22 «23 «24 
«31 «32 «33 «34 
\ «41 «42 «43 «44 / 



(2.13) 



acting on nearest neighbour sites j and j + e*-"^ . The sum runs over the whole lattice and 
we take periodic boundary conditions in all space directions, i.e., ^(a) j{d)) acts 

on sites ■ ■ ■ ,L^"'\ . . . , j^'^-') and ■ ■ ■ , . . . The normalization sets the 

time scale and is for the purposes of this paper of no particular interest. The diagonal 
elements a^k of satisfy 

4 

cifefc = - X! (2-14) 

k' = l 

which is imposed by conservation of probability. This implies 

(s|ui"^=0 Vj,a (2.15) 

In order to keep the interpretation of H as defining a stochastic process, (|1|) has to be 

supplemented by the condition akk' ^ for the off diagonal matrix elements k ^ k' . 

The processes decribed by H are reactions changing the configurations on two nearest 

neighbour sites. A configuration {nj,njj^g(a)} changes into configurations {n'jjn'j^^^^} 

with rates akk' as follows: 

{0, 0} a2i{0, 1} + asijl, 0} + a4i{l, 1} (birth/pair creation) 

{0, 1} ai2{0, 0} + a32{l, 0} + a42{l, 1} (death/diffusion/decoagulation) 

(2.16) 

{1, 0} —>■ aiajO, 0} + a23{0, 1} + 043(1, 1} (death/diffusion/decoagulation) 

{1, 1} ai4{0, 0} + a24{0, 1} + 034(1, 0} (pair annihilation/coagulation) . 

These processes take place with equal rates everywhere in the lattice. 

This generalized nearest neighbour exclusion process includes many well-known pro- 
cesses such as the asymmetric exclusion process |l^ (with hopping rates 023,032 7^ 0, 
all other rates 0), the voter model ||l^ (with death rates and decoagulation rates 
012 = oi3 = 042 = 043 7^ 0) or Glauber dynamics |l^] (023 + 032 = 014 + 041 / 0). 
Altogether there are independent 12 parameters, one of which is trivial as one may always 
change the normalization ^ without changing the physical properties of the system. We 
shall set ^ = 1 throughout the paper. 

3 Equal-time correlation functions 

The equal-time fe-point correlation function satisfies the equation 

^( • • • n^, ) = - ^ J2 ( • • • n^.-ui"^ ) (3.1) 

a=l j(zc(<^) 
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where, owing to the property (2.15) of the two-site reaction matrix the sum over j does 

not run over the whole lattice, but only over the union C^"^ of the set of sites {xi, . . . , Xk} 
with the set of their nearest neighbours {xi — e^^\...,xl. — e*-"^}. It is important to 
realize that since each Uj^^ acts non-trivially only on two sites, the r.h.s. of ( |3.l| ) involves 

only (k — 2)-point functions, (k — l)-point functions, A:-point functions and {k + l)-point 
functions. This can be seen as follows: Suppose one of the Xi G {xi, . . . ,Xk} (say Xk) is 
equal to j + e^""^ . Using ( |2.10| ) one finds 



(3.2) 



with 



(3.3) 



M = a2l + a4i Bi = 023 + 043 - ^21 - 041 

Cl = ai2 + 032 + 021 + 041 Di = Ci - 023 " ^43 " OU " 034 • 

A similar result arises if one of the xi (again, without loss of generality Xk) is equal to 



with 



( • • • ("-xfc 4"^ ) ) = ^2 ( • • • n^^.i ) 
M = 031 + 041 B2 = 032 + 042 - 031 - 041 

C2 = Oi3 + O23 + 031 + 041 D2 = C2- 032 - 042 - Oi4 - 024 



(3.4) 



(3.5) 

The r.h.s. of ( |3.2| ) and ( |3.4| ) consist only of (A; — l)-point functions, /c-point functions and 
{k + 1) -point functions. 

If two of the Xi are nearest neighbours in the lattice, e.g. Xk-i = Xk — e^"^ = j, then 
the action of Uj''' yields 

( • • • ^x,_2(nx,_i%^_^_^^-<„)4t^_ J ) = A^{ns^ ■ ■ ■ n^^_2 ) 



with 



A3 = 041 B3 = 042 — 041 

-^3 = O43 — O41 C3 = Oi4 + O24 + O34 + O42 + O43 — O41 . 



(3.6) 



(3.7) 
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The r.h.s. of ( |3.6D consists only of {k — 2)-point functions, {k — l)-point functions and 
/c-point functions. 

li Di or D2 7^ the time derivative ( ^.iD of the A;-point correlation function gives 
rise to a set of M coupled differential-difference equations involving all fc-point correlators.]^ 
Solutions to such a set of equations have been found in some special cases where subsets 



of these equations decouple |21|, but there is no general solution. On the other hand, if 
Di = D2 = 0, i.e., for 

034 = 021 + 041 + ai2 + 032 - a23 - 043 - ai4 

(3.8) 

124 = asi + 041 + ai3 + 023 — ^32 — 042 — 014 , 

the problem simplifies considerably as one has a closed system of only k equations. In this 
case Eq. ( |3.lD may be regarded as an inhomogeneous, linear differential-difference equation 
for the A;-point function with (/c — l)-point and (/c — 2)-point correlators as inhomogeneities. 

Note that one may also study correlation functions of the operators h = n — a with 
an arbitrary constant a. One obtains again a closed system of k equations for the A;-point 
correlation function if Di = D2 = 0, but with new constants 



and 



Ai = Ai + a{Bi-Ci) , Ai = Ai+a{Bi-Ci) (3.9) 

ig = A^ + a{B^ + D^)-a^C^ 

= Bs - a{Cs + B2 - Ci) (3.10) 
Ds = 53-0(^3 + ^1-^2) 



in Eqs. (p^), ( |3.4[) , (|3.6|). i?i_2 and Ci^2 do not change. In particular, the inhomogeneity 

with 

(3.11) 



arising from Ai,A2 7^ in the one-point function can be removed by taking a = p with 

2a4i + 021 + 031 Ai + A2 



2a4i + 021 + 031 + 2ai4 + 024 + 034 Ci + C2 — Bi — B2 

With this choice one has Ai + A2 = 0. The differential-difference equation for higher 
order correlation functions have then inhomogeneous terms proportional to ^3 (coupling 
to {k — 2)-point functions) and S3 = S3 + Ai, D'^ = + A2 (coupling to {k — l)-point 
functions. We conclude: 



Eq. I[3. l\) becomes a closed, inhomogeneous, linear differential- difference equation in one 
(continous) time coordinate and d ■ k (discrete) space coordinates on a 10-parameter sub- 
manifold defined by Eqs. ( \3.S{ ) of the 12-parameter model. 

One may add the remark that this differential-difference equation becomes homogeneous 
(i.e., contains no {k — l)-point and (k — 2)-point correlation functions) on a 7-parameter 
submanifold defined by A3 = B'^ = D'^ = 0. From the derivation presented above it 
is obvious that this result is easy to generalize to other lattices and interactions. With 
Di = D2 = Q Eq. ( |3.1[) is a closed set of k equations independent of the dimensionality 
of the system or of the kind of lattice on which the model is defined. Furthermore, the 
two sites on which the reaction matrix u acts nontrivially are not even required to be 

^In A^-particle systems with particle number conservation the hierachy breaks off at A; = A'' < M. 
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nearest neighbours. The result remains true for arbitrary long-range interactions with 
reaction matrices u^^^ where x and y are any two points on the lattice. Finally, it is 
also not necessary to keep the reaction rates akk' space-independent as long as (3.8) is 
satisfied for each reaction matrix u^g. For the decoupling from lower order correlators 
a stronger condition is necessary for this general case. Besides Di = D2 = one needs 
Ai = A2 = A3 = B3 = D3 = 0, i.e., one is left with a five-parameter space only. 

We demonstrate this result for the one-dimensional case. In one dimension, Eq. ( |3.lD 
for the one-point function becomes 

d 

g^irix) = -{rixiux-i +Ux)) (3.12) 

where we have set ^ = 1 and dropped the coordinate index a in This is easy to 

calculate and one finds 



^{rix) = A1+A2 



+Bi{nx-i)- {Ci + C2){nx) + B2{nx+i 
+Di{nx-inx) + D2{nx rix+i) 



(3.13) 



where the constants with index 1 and index 2 arise from the action of u^-i and Ux re- 
spectively. One sees that if Di = D2 = ( p^ ) b ecomes an inhomogeneous, linear 
differential-difference equation. Introducing n with a = p as defined in ( |3.1l| ) leads to the 
homogeneous equation 

^^{hx) = Bi{fix_i) + B2{hx+i) - {Ci + C2){nx) . (3.14) 

For the two-point function one obtains 
d 

■^{rixny) = -{nx{ux-i +Ux)ny) - {nxny{uy-i +Uy)) 

= {Ai+A2){{nx) + {ny)) ^^-^^^ 
Bi{{nx-iny) + {nxUy^i)) + B2{{nx+iny) + (n^%+i)) 
-2{Ci + C2){nxny) 

if X and y are not nearest neighbours and 
d 

— {UxUx+l) = -{UxUx+liUx-l +Ux +Ux+l)) 

= A3 + {A2 + D3){nx) + {Ai + B3){nx+i) ^g^^g^ 
+Bi{nx-inx+i ) + B2{ 
-{Ci + C2 + C3){nxnx+i) . 

for the nearest neighbour correlator. Similar equations are obtained for correlators involv- 
ing n. 
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For Di = D2 = and the special choice aki = = (no birth, pair creation 
and decoagulation) one has Ai = = D3 = and ( p.l6| ) simphfies to the completely 
decoupled, homogeneous equation 

d 

— (nj.na;+i) = Si(nx._in^+i ) + -B2(na;n^.+2 ) (3 17) 

-(Ci + C2 + C3)(n^n^.+i) . 

In this case, time derivatives of higher fc-point correlation functions also decouple com- 
pletely, i.e., involve only /c-point correlation functions. 



4 Dual processes 

The fact that the time-derivative of the A;-point correlation functions may give rise to a 
closed set of equations is reminiscent of the duality relations e.g. for the symmetric exclu- 



sion process [15, IC] where self-duality is indeed just an expression of this fact. One may 
therefore ask whether the closure of the equations for the /c-point functions is equivalent 
to the existence of some dual process. 

Before we discuss this question we would like to remind the reader of the meaning of 
self-duality (of the symmetric exclusion process) in the operator language used in this pa- 
per. Let us assume that initially N particles are located on a set of sites An = {yi, . . . , yXr} 
represented by a vector \Aj\f) = \yi, . . . ,yp^). We want to compute the probability 
( n^^ • • • n^j. )am of finding (any) k particles on sites = {xi, . . . ,Xk}, at time t. The 
duality relations state ||l5|, |9| 



AT 



In this expression the sum runs over all sets 5^ = {x'l, . . . , x^} which are contained in the 
set A]^, i.e., the fc-point correlation function (n^^ • • •n^^)Aj^ of the A''-particle system is 
given by sums of fe-particle correlation functions (we assume A; < iV). Using (l2lc| ), ( ^ll|) 



and the fact that H for the symmetric exclusion process is symmetric and SU (2) invariant 
(i.e., commutes with S^), the duality relations ( ^T| ) can be derived as follows j^: 



{nn^-'-n^jAN = {s\ns-,■■■n^^^e ^^\An) 



Y.i^N\e^~\n){xi,...,Xk\e-"'\n) (4.2) 

n 

Y: {x[,...,x',\e-"^'\Bk) 
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where \Bk) = . . . Because of particle number conservation we have substituted 

{x'l, . . . ^x'j^ \ in the last line by ( s \nff^ ■ ■ ■ n^i . In other words, the averaging is performed 
over all fe-particle states such that the sets of occupied sites are contained in the set 
A]\j of initially occupied sites. These sets arise from the matrix element {A]\[\e^ \n) 
together with particle number conservation. 
Generally we define a duality relation by 



{s\Qe-^'\A) = {s\Q'e-"'\A' 



(4.3) 



where Q and Q' are some functions of the projection operators and | A ) and | A' ) are 
initial states. The dual process H is obtained by taking the transposed matrix of the 
time evolution operator and performing some suitably chosen similarity transformation V 
such that 

H = VH'^V-^ (4.4) 



are 



indeed defines a stochastic process. The observable Q' and the initial condition | A' ] 
then given by 

{s\Q' = {A\V-^ , \A') = VQ\s) (4.5) 

where | ...) = ((.. . |)'^. By taking V = exp {—S~) and Q = n^^ ■ ■ ■ n^^ one recovers ( [4.1[) 
with H = H. 

In this kind of duality the points {xi, . . . ,Xk} occuring in the correlator are mapped 
to an initial state with particles occupying sites {xi, . . . ,Xfc}. The dual process may be 
considered as a process describing the time evolution of particles on these points. H = H 
means that the symmetric exclusion process is self-dual. 

After this reminder we are in a position to formulate the problem more specificly: 
We have seen in the previous section that as in the symmetric exclusion process the 
time derivative of the fe-point function of the generalized model does not involve higher 
correlators. The transformation V = exp (—S~) relates the A:-point correlator to a fc 



particle initial state. We therefore ask the question whether the dual operator H (iA) with 
this particular V defines a stochastic process on the 10-parameter manifold Di = D2 = 0. 
In order to answer this question one has to check whether the transformed dual rates a^i 
satisfy (|2.14| ) (guaranteeing conservation of probability) and the condition of positivity 
fflfei > for 7^ L 

A short calculation gives for the dual matrices 



-Xa) 



( 







A2 


^3 \ 


Ci 


B2 




Bi 


-C2 










-C3 ) 



(4.6) 



0» 



with the quantities Ai, Bi, Ci defined by ( |3.3D , ( ^ ) and ( |3.7D respectively. Using positivity 
and conservation of probability the conditions for the existence of a dual process can be 
read off, namely Ai,Bi,Dz > and Ai + Bi-Ci = A2 + B2-C2 = AS + B3 + D3-C3 = 0. 
Note that the condition C3 = 042+043 — 041+014+024+034 = A3+B3+D3 = 042+043 — 041 
implies 0^4 = OVA; because of the positivity of the original rates. 
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The dual process H has only hopping terms Bi and B2 and various non-vanishing 



annihilation rates, but no particles are created. Thus the duality relations (O) with 
Q = n^i--- "-Xfc read 



+ •••+ E {yi,...,yk\e-'''\Bk) (4.7) 

s7i,...,s7fcGAjv 

k 

= E E {^\Q'ip)\Bk) 

P=OB'pCAn 

where Q'{p) = Ppnyi ■ ■ ■ n^- and Pp is the projector on p-particle states arising from the 
matrix element {A^le^ | n). The averaging extends therefore over all p-particle states 
with < p < k such that the sets B' of occupied sites are contained in the set of 
initially occupied sites. Note that ( [4.7|) holds for any choice of the parameters Ai, Bi, D^, 
irrespective of whether H defines a stochastic process or not. 

One may define other dual processes involving other operators Q. The main require- 
ment for the kind of dual processes we are interested in is that the set of sites B^ defined 
by the product of projectors rig- translates into an initial configuration of occupied sites. 
This feature determined the transformation V. The results of the preceding section in- 
dicate that it might be interesting to study correlation functions of the shifted density 
projectors We introduce the local operator 

(1 + a a \ 
(4.8) 
—a 1 — a I ^ 
' 3 

and 

W = ^w.. (4.9) 

J 

It is easy to check that {s\w^ = {s \ and ( s |nj = {s\ (rej — a)wj. Therefore 

( s \ng^ ■■■hg^ = {s \ns^ ■ ■ ■ n^^W (4.10) 

and 

{n£^---ns^)A = {s\n^^---ng^(r^'^\A') (4.11) 

with H' = WHW^^ and the transformed initial state \ A' ) = W\ A). Now one may apply 
(O) to diHl ). Avera ging is now performed over the weighted set of states n with total 



particle number p < k and with weights fA{yi,---,yk) = ( ^' |e'^ \yi, . . . ,yk)'- 

{ns^- ■ -n^k) A= E fA{yi,---,yk){yi,---,yk\e~"*\Bk) . (4.12) 



0<p<k 
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One obtains the dual time evolution operator H 
transformed dual reaction matrices 



V{WHW-^YV-^ with the doubly 



/ 


Ai 


A2 


is \ 





-Ci 


B2 


S3 





Bi 


-C2 


^3 


V 








-C3 ) 



(4.13) 



0» 



where the quantities defined in Eqs. ( |3.9D , ( |3.10| ) are used. Positivity and conservation of 
probability yield again the conditions on the existence of the dual process defined by H. 
We conclude that 

the closure of the differential- difference equations for the k-point density (or shifted den- 
sity) correlation functions does not, in general, imply the existence of a dual stochastic 
process as discussed here. Additional constraints on the original reaction rates arise in 
order to conserve probability and positivity for the dual process. 

Other duality relations and dual processes may be obtained by considering other correla- 



tion functions |12] or transformations to other initial states. 



5 The Heisenberg Hamiltonian 



In the previous section we have shown that on a 10-dimensional submanifold of the 
12-parameter problem all equal-time correlation functions can be calculated by solving 
(in)homogeneous, linear differential-difference equations. From the solutions to these equa- 



tions one may obtain the spectrum of H (2.12) by looking for the poles of the Laplace 
transform of the correlation function. From the one-point function (3.1) or (6.1) below 
one finds a series of eigenvalues 



a 

E{k) = [Bie'^-^ + B2e-'^^ -C1-C2) < ka < 2tt 



(5.1) 



a=l 



which may be interpreted as non-relativistic, free single-particle excitations. The full 



spectrum would be obtained from the solution to all correlation functions. From (5.1) we 
find that the system is partially integrable in any number of space dimensions. Partial 
integrability is known to occur also in a different, 7 parameter, subspace of the model in 



one dimension |21|, but there the energies have a more complicated structure. Note that in 
the discussion in the previous section the positivity of the constants ati was only necessary 
for the interpretation of H as generating as stochastic process. The partial integrability 
of H is ensured by the constraints ( |3.8D alone. 

In order to get some insight into the physical origin of the one-particle excitations and 
of the structure of the equations for the higher order correlators we study the relation- 
ship of the stochastic Hamiltonian ( 2.12)) to the Heisenberg quantum Hamiltonian defined 
below. This is motivated by the symmetric exclusion process in which case (|1|) is the 
Hamiltonian of the isotropic Heisenberg ferromagnet. 
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In this section we show that 



the spectrum of H ( 2.1^ with the constraint jjjS.S^ is identical to the spectrum of the Hamil- 
tonian of an anisotropic spin- 1/2 Heisenberg quantum Hamiltonian Hxxz in o- magnetic 
field. The behaviour of the k-point correlation functions is determined by the excitations 
of the p-magnon sector where < p < k. 



The ferromagnetic Heisenberg quantum Hamiltonian Hxxz is of the same form as 
( 2.12| ), but with matrices 



h 



(a) 



2-v V : 



/ 











\ 







1 










1 









Vo 








/l44 


/ 



(5.2) 



0» 

where have chosen the normahzation 7=1 and the constant c = — (61 + 62 + ^) such 
that /i23 = h^2 = 1, hii = 0, /i22 = -(/32 + A), /133 = -(/3i + A), /144 = - /32, 
and /i23 = /132 = 1- This Hamiltonian has a continous U{\) symmetry generated by 
5^ = Hence Hxxz sphts into sectors with fixed z-component of the total spin. 

By identifying spin up at site j with a vacancy and spin down with a particle, this U{1) 
symmetry amounts to particle number conservation. The sector with 5^ = M/2 — N 
corresponds to the A/'-particle sector. 

In order to make contact with the reaction matrices Uj^^ we study the dual matrices 

u'f^ ( [4. 61) which where obtained from the original reaction matrices by the similarity trans- 
formation V and transposition (|4.4| ). First we renormalize -ff by a factor ^ = ^BiB2 and 
perform another similarity transformation H = {^H^^^)/^ where $ = exp " 

and ff = r] Ea=i e^''^ with r/ = y/W/B^- For ii = i2 = ^3 = i^3 = -D3 = this 
manipulation yields matrices 

/ \ 

U22 1 
1 n33 

\ U44 y 



(5.3) 



0» 



with Till = 0, n22 = Ci/S,, = C2/C, ^44 = U32 = ^23 = 1- These matrices are iden- 
tical to the Heisenberg matrices ( |5.2| ) for an appropriate choice of Pi, (32 and A. Therefore 
on the five-parameter submanifold of the general model on which one obtains completely 
decoupled equations for the fc-point correlation functions one finds H = Hxxz^ The 

^The transformation $ relates the symmetric model /123 ~ h^2 to the asymmetric model ^23 7^ '^32, see 
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duality transformation turns the state ( s \n^^ ■ ■ ■ n^^ ( [4.10| ) appearing in the /c-point cor- 



relation function (4.11) into a A;-particle initial state. Since Hxxz (and the transformed 
stochastic Hamiltonian H) conserve particle number, one finds that the correlator is given 
by the /c-magnon excitations of the Heisenberg quantum Hamiltonian. 

In order to understand the more general 10 parameter manifold Di = D2 = we split 
the transformed dual reaction matrices into two parts u 



(a) 



h 



(a) 



/I'l"^ with 



h' 



(a) 






il 


A2 


is 











^3 











D3 















(5.4) 



0» 



Thus one may write H = Hxxz + H'. The crucial point is that Hxxz conserves particle 
number and may therefore be block diagonalized in blocks with fixed N . On the other 
hand, H' connects a block with particle number with blocks with particle numbers A'^— 1 
and N — 2, but not with any N' > N. The whole matrix H has therefore a triangular 
structure with block matrices labelled by N arising from Hxxz on the diagonal and matrix 
elements resulting from H' on the upper off-diagonal. The characteristic polynomial of H 
does not depend on these off-diagonal entries and therefore the characteristic polynomials 
of the stochastic Hamiltonians H, H (|2.13| ),( ^!8|) and the Heisenberg Hamiltonian Hxxz 
(5.2) are identical.^ The full dynamics of the A;-point correlation function are determined 
not only by the eigenvalues of the time evolution operator, but also by the eigenstates. 
Since H' annihilates only, eigenvalues and eigenstates with p < k particles contribute 
to the dynamics. This explains the partial integrability of the model: To the one-point 
function k = 1 only the one-magnon sector contributes. This is indeed a non-relativistic 
free particle. It is interesting to note that in one dimension, Hxxz is completely integrable. 
In this case the eigenvalues may be found from the Bethe ansatz. 



6 The average density 

In this section we give an application of the results of Sec. 3. The simple form of Eq. ( |3.1D 
for the shifted average particle density on a hypercubic lattice, 



Ql{f^x) = Y. %--et^) > + B2{fl^+^a) ) - (Cl + C2){h^ )) , (6.1) 

a=l 

allows for an explicit integration and thus the extraction of the critical and non-critical 
behaviour of the system. We define 

C = Ci + C2, = 2^/3^2, '^ = ^ln(^|^) (6.2) 



^ [2^, Q for the one-dimensional case. It induces non-periodic, twisted boundary conditions, i.e., the 
constants 1123 and 1132 in the boundary matrices are different from those in the bulk matrices u^f \ In so 

far Hxxz and H agree in general only up to boundary terms. 

^ The same argument was used in Ref. , but on another submanifold of the parameter space. 
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and the vector £ = £ Yl'a=i ^""^ study first the infinite system. In d dimensions the 
solution p{x,t) = {hg) to (|6.lD is in terms of modified Bessel functions In{x) given by 



p{x, t) = Y: ag e-'<^^+^<'-y^ n 4_^^ {Dt) (6.3) 
y a=i 

with initial condition 0) = ag. With = adg^^ and for large times t and large 
separations = {x — y)"^ with /t fixed this may be written 

p(r, t) = ae-'^^*(27rZ?t)-'^/2g-i^^^ (g_4) 

and the physical interpretation of the constants becomes apparent. D plays the role of a 
diffusion constant, while i5 is a driving field leading to an average drift velocity v = D£. 
From this we find that the system satisfies the Einstein relation 



'd£^ 



I =D . (6.5) 

£:(")=o 



The constant 

P = C-{l + ^-^)D (6.6) 

is a decay constant. For /i = the system is critical with dynamical exponent z = 2. 

For a finite-size scaling study of the density on a hypercubic lattice we consider a one- 
dimensional system with L sites. The solution to ( |6.1| ) with periodic boundary condition 
and initial condition /j(x, 0) = adx^y (with 1 < x,y < L and r = x — y) is 

Kr, t) = ^ exp I — ^ [2C - {Bi + ^3) cos — - i{Bi - B2) sin — ]t j (6.7) 

Introducing the scaling variables u = 'i:{r + [Bi — B2)t)/L and r = 27r(i?i + B2)t/L? and 
taking the limit L ^ 00 this becomes 

p(tx,r) = ^e-^'*03(n|ir) (6.8) 

with the decay constant jj,' = C — Bi — B2 and the Jacobi theta-function 

00 

^3(n|ir) = e-™'^+2^'^" . (6.9) 



Eqs. (6/7) and (|6.8D describe a density distribution of width A = i?i + ^2 = i^cosh£^ 
with its center aX x = y — vt where v = Bi — B2 = D sinh^". The Einstein relation ( |6.5D 
therefore hold also in the finite system. Note that when taking the scaling limit we have 
implicitly assumed that £ = e/L vanishes proportional to L"^ because otherwise u would 
diverge. This implies A ~ D, ~ D£ and fx' ^ fi. For /i' = and large values of r, i.e. for 
times larger than L^, the density in the moving frame of reference decays exponentially to 
its constant value poo = cl/L with decay constant 27r^-D. For small values of r the density 
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decays with a power law behaviour, /j(0, r) = aL ^ (27r -Dr) ^Z^. This may be shown using 
the Poisson resummation formula, but is already clear from (|6.4| ). 

Eq. Q may be solved for other boundary conditions such as open boundary con- 
ditions with injection and absorption of particles [25|. We do not want to discuss this 
here. 

If j4i, ^2 7^ the constants Bi and B2 can be negative. On bipartite lattices this leads 
to an alternating positive and negative (relative) density p{x,t). For other lattices the 
situation is more complex. 



7 Conclusions 

We have studied a general reaction diffusion process of hard-core particles with two-site 
interactions on a lattice in d dimensions. There are twelve parameters ( p.l6|) for the 



various reaction and diffusion rates. On a 10 parameter submanifold defined by (|3.8|) 



of the parameter space the differential-difference equations satisfied by the /c-point 
density correlation functions are inhomogeneous linear equations involving no higher order 
correlators. We have mainly considered a d-dimensional hypercubic lattice with nearest 
neighbour interaction and space-independent rates, but this result holds for arbitrary 
lattices with arbitrary two-site interactions as long as (^^) is satisfied for the interaction 
between each pair of sites. On a d-dimensional hypercubic lattice with periodic boundary 
conditions Eqs. ( ^.l]) decouple completely on a 7 parameter submanifold. 

It is perhaps worthwhile pointing out that throughout Sections 3 - 5 we have assumed 
that initially particles are located on a set of sites {xi, . . . ,xn}. All calculations can 
be easily repeated for an arbitrary, time-dependent initial distribution. In this way one 
can get similar results for two-time correlation functions. 

We have shown that in general the decoupling from higher order correlators does not 
imply the existence of a dual process for the time evolution of a /c-particle initial state as 
one has in the special case of the symmetric exclusion process. This remains true only 
under further assumptions on the reaction rates arising from conservation of probability 
and positivity of rates in the dual process. The dual process to the general 10 parameter 
model (if it exists) is a process involving only diffusion and annihilation of particles (death, 
decoagulation and pair annihilation), but no creation. 



The Hamiltonian (2.13) defining the stochastic process is partially integrable and has 
the same spectrum as a spin-1/2 Heisenberg quantum Hamiltonian Hxxz in a magnetic 
field and with twisted boundary conditions. This explains the occurence of the non- 
relativistic one-particle excitations appearing in the time evolution of the density profile 
in terms of one-magnon excitations of the Heisenberg model. The dynamics of the k- 
point correlators are given by p magnon states with p < k. In one dimension, Hxxz is 
completely integrable and its spectrum can be found from the Bethe ansatz. However, it 
is interesting to note that even though for Di = D2 = the spectra of H and Hxxz stice 
identical, the matrices Uj do not, in general, satisfy the Hecke algebra relations UjUj±iUj — 
Ui = Uj±iUjUj±i — Uj±i, [ui,Uj] = for |i — j| > 2 and ii| = Xuj [13, U]. Through 
Baxterization |22| this would imply the integrability of the model as in the case of the 
normal Heisenberg Hamiltonian ( |5.2| ). This observation hints at a more general algebraic 
structure beyond the usual conditions for integrability in one dimension. 
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In Sec. 6 we have used our results for the exact calculation of the time evolution of 
the density from an arbitrary initial density and analyzed its finite-size scaling behaviour 
in the scaling regime t ~ (L is the size of the system). An initially sharp peak in 
the distribution widens diffusively and moves with a constant average velocity. It turns 
out that the model satisfies the Einstein relation ( |6.5D relating the drift velocity and the 
diffusion constant. Depending on the various reaction rates there is a non-critical region 
with an additional exponential decay of the amplitude. Thus the time evolution of the 
density profile on a d-dimensional hypercubic lattice with periodic boundary conditions 
depends only on four combinations of the 10 parameters. For other boundary conditions 
or lattices more parameters enter and it would be interesting to study the corresponding 
lattice effects. The lattice diffusion constant may be or even negative. The latter case 
corresponds to the development of an alternating structure of the average density on 
bipartite lattices before reaching the constant stationary density. Another open question 
is the behaviour of the density-density correlation function in the general 10 parameter 
model. This quantity may be non-trivial even in the steady state as only under strong 
restrictions on the rates the steady state is a product measure. 
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